CLNS 08/2031 
CLEO 08-14 

Observation of Y(2S) — ► r]X(lS) and search for related transitions 

Q. He, 1 J. Insler, 1 H. Muramatsu, 1 C. S. Park, 1 E. H. Thorndike, 1 F. Yang, 1 
M. Artuso, 2 S. Blusk, 2 S. Khalil, 2 J. Li, 2 R. Mountain, 2 S. Nisar, 2 K. Randrianarivony, 2 
N. Sultana, 2 T. Skwarnicki, 2 S. Stone, 2 J. C. Wang, 2 L. M. Zhang, 2 G. Bonvicini, 3 
D. Cinabro, 3 M. Dubrovin, 3 A. Lincoln, 3 P. Naik, 4 J. Rademacker, 4 D. M. Asner, 5 
K. W. Edwards, 5 J. Reed, 5 R. A. Briere, 6 T. Ferguson, 6 G. Tatishvili, 6 H. Vogel, 6 
M. E. Watkins, 6 J. L. Rosner, 7 J. P. Alexander, 8 D. G. Cassel, 8 J. E. Duboscq*, 8 
R. Ehrlich, 8 L. Fields, 8 R. S. Galik, 8 L. Gibbons, 8 R. Gray, 8 S. W. Gray, 8 
D. L. Hartill, 8 B. K. Heltsley, 8 D. Hertz, 8 J. M. Hunt, 8 J. Kandaswamy, 8 
D. L. Kreinick, 8 V. E. Kuznetsov, 8 J. Ledoux, 8 H. Mahlke-Kriiger, 8 D. Mohapatra, 8 
P. U. E. Onyisi, 8 J. R. Patterson, 8 D. Peterson, 8 D. Riley, 8 A. Ryd, 8 A. J. SadofT, 8 
X. Shi, 8 S. Stroiney, 8 W. M. Sun, 8 T. Wilksen, 8 S. B. Athar, 9 R. Patel, 9 J. Yelton, 9 
| 7 j ! P. Rubin, 10 B. I. Eisenstein, 11 I. Karliner, 11 S. Mehrabyan, 11 N. Lowrey, 11 M. Selen, 11 

E. J. White, 11 J. Wiss, 11 R. E. Mitchell, 12 M. R. Shepherd, 12 D. Besson, 13 T. K. Pedlar, 14 
J. V. Xavier, 14 D. Cronin-Hennessy, 15 K. Y. Gao, 15 J. Hietala, 15 Y. Kubota, 15 
T. Klein, 15 B. W. Lang, 15 R. Poling, 15 A. W. Scott, 15 P. Zweber, 15 S. Dobbs, 16 
V : Z. Metreveli, 16 K. K. Seth, 16 A. Tomaradze, 16 J. Libby, 17 L. Martin, 17 A. Powell, 17 

G. Wilkinson, 17 K. M. Ecklund, 18 W. Love, 19 V. Savinov, 19 H. Mendez, 20 J. Y. Ge, 21 
D. H. Miller, 21 I. P. J. Shipsey, 21 B. Xin, 21 G. S. Adams, 22 M. Anderson, 22 
J. P. Cummings, 22 I. Danko, 22 D. Hu, 22 B. Moziak, 22 and J. Napolitano 22 

(CLEO Collaboration) 

1 University of Rochester, Rochester, New York 14627, USA 
• 2 Syracuse University, Syracuse, New York 13244, USA 

\£2 '. 3 Wayne State University, Detroit, Michigan 48202, USA 

4 University of Bristol, Bristol BS8 1TL, UK 
5 Carleton University, Ottawa, Ontario, Canada K1S 5B6 
6 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
7 Enrico Fermi Institute, University of Chicago, Chicago, Illinois 60637, USA 
^ \ 8 Cornell University, Ithaca, New York 14853, USA 

9 University of Florida, Gainesville, Florida 32611, USA 
10 George Mason University, Fairfax, Virginia 22030, USA 
11 University of Illinois, Urbana- Champaign, Illinois 61801, USA 
12 Indiana University, Bloomington, Indiana 47405, USA 
13 University of Kansas, Lawrence, Kansas 66045, USA 
^Luther College, Decorah, Iowa 52101, USA 
15 University of Minnesota, Minneapolis, Minnesota 55455, USA 
16 Northwestern University, Evanston, Illinois 60208, USA 
17 University of Oxford, Oxford 0X1 3RH, UK 
18 State University of New York at Buffalo, Buffalo, New York 14260, USA 
19 University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 
20 University of Puerto Rico, Mayaguez, Puerto Rico 00681 



> 



oo 
o 



* Deceased 



1 



21 Purdue University, West Lafayette, Indiana 47907, USA 
22 Rensselaer Polytechnic Institute, Troy, New York 12180, USA 

(Dated: June 17, 2008) 

Abstract 

We report the first observation of T(2S) — > wY(lS), with branching fraction B = (2.1+Qg (stat.)± 
0.3 (syst.)) x 10~ 4 and statistical significance 5.3a. Data were acquired with the CLEO III detector 
at the CESR e + e _ symmetric collider. This is the first process observed involving a 6-quark spin 
flip. For related transitions, 90% confidence limits in units of 10~ 4 are B[T(2S) — > 7T°T(1S')] < 1.8, 
B[T(3S) -► nT(lS)] < 1.8, B[T(3S) -► tt°T(LS)] < 0.7, and B[T(3S) -► ^°T(2,S)] < 5.1. 
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In order to produce a pseudoscalar meson r\ or 7r° in T(nS) — > (rj/7r°)T(mS) transitions, 
the bb pair must emit either two Ml (chromomagnetic dipole) gluons or an El (chromoelectric 
dipole) and an M2 (chromomagnetic quadrupole) gluon [l|, 0,0], involving the flip of a heavy 
quark's spin. In this Letter we present the first observation of T(2S) — > r]T(lS), and a search 
for similar tt° or 77 transitions from the T(2S) and T(3S). A spin-flip of a 6-quark can shed 
light on its chromomagnetic moment, expected to scale as 1/m.f,. Electromagnetic transitions 
involving a 6-quark spin-flip should also have amplitudes scaling as l/m^. They have not 
previously been observed. 

The decay ip{2S) — > T]J/ip was observed in the early days of charmonium spectroscopy 
0]. Its branching fraction is B[ijj(2S) — > i]J/ip] = (3.13 ± 0.08)% while only an upper 
limit B < 2 x 10~ 3 is known for the corresponding T(2S) — > 1]T(1S) process 0. The 
upper limit for T(35) -> r]T(lS) is B < 2.2 x 10~ 3 Q. The quark spin-flip involved in 
T(nS) — > (77 / 7T°) T(mS) transitions (we consider 3 > n > m > 1) and the P-wave nature of 
the final state imply that rates should scale from charmonium as T oc (p*) 3 / m g [HQ]; where 
p* is the three-momentum of the 7] or tt° in the T(nS) center-of-mass system and Q = c,b 
is the heavy quark. Hence one expects 

r[T(2S,3S)^T(lg)] _ 

Ity(2S) - ijJ/y>] -^- uu ^ u - UUi ^' W 

leading to £[T(2S,3S) -> r/T(15)] ~ (8.0,6.5) x 10^ 4 . Direct calculation in a potential 
model [2| yields (6.9, 5.4) x 10 -4 for these branching fractions. All predictions involve a 
perturbative calculation of gluon-pair emission followed by a nonperturbative estimate of 
materialization of the gluon pair into an 77. Uncertainties associated with this estimate are 
difficult to quantify. 

Similar predictions can be made for tt° transitions under the assumption that they are 
due to an isospin-zero admixture in the 7r°. The isospin-forbidden decay ip(2S) — > tt J/i/j 
has been seen jBJ with a branching fraction of (1.26 ± 0.13) x 10~ 3 which is (4.03 ± 0.43)% 
of that for ip(2S) — > rjJ/ip. Using values of p* appropriate to each process and assuming the 
same isospin-zero admixture in 7r° governs the transitions T(nS) — > ir°Y(mS), one obtains 
the scaling predictions 

B[T(2S,3S)^n°T(lS)] 

B[T(2S,3S)^ V T(1S)] = (16 ± 2 ' °- 42 ± °-° 4)% • (2) 

There is no prediction at present for the kinematically-allowed decay T(35) — > tt°T(2S). 

The data in the present analysis were collected in e + e~ collisions at the Cornell Electron 
Storage Ring (CESR), at center-of-mass energies at and about 30 MeV below the T(2S, 3S) 
resonances. Integrated luminosities at these resonances were (1.3,1.4) fb _1 , amounting to 
(9.32±0.14, 5.88±0.10) million decays of T(2S, 3S), as in the analysis of Ref. [8j. Events were 
recorded in the CLEO III detector, equipped with an electromagnetic calorimeter consisting 
of 7784 CsI(Tl) crystals and covering 93% of solid angle, initially installed in the CLEO II 
[s| detector configuration. The energy resolution of the crystal calorimeter is 5% (2.2%) for 
0.1 (1) GeV photons. The CLEO III tracking system [10] consists of a silicon strip detector 
and a large drift chamber, achieving a charged particle momentum resolution of 0.35% (1%) 
at 1 (5) GeV/c in a 1.5 T axial magnetic field. 

We look for candidate events of the form e + e~ — > T(nS) — > (rj / 'tt°)T \mS) with T(mS) — > 
£ + £~, where I = e, /i. Candidates for are identified by picking the two highest-momentum 
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tracks in an event and demanding them to be of opposite sign. We explore separate e + e~ and 
samples in T(mS) decays by defining electron candidates to have a high ratio of energy 
E observed in the calorimeter to momentum p measured in the tracking system, i.e., E/p > 
0.75, and muon candidates to have E/p < 0.20. We choose lepton candidates from tracks 
satisfying | cos#| < 0.83, where 9 is the angle with respect to the positron beam direction, 
to avoid a region of less uniform acceptance at larger | cos#|. With these criteria we achieve 
a very clean separation of electron and muon candidates. In order to suppress contributions 
from Bhabha scattering, we demand for events with (r], tt°) — > 77 that e + candidates satisfy 
cos# e + < 0.5. This greatly suppresses Bhabha scattering background while keeping 93% 
of the signal. Once leptons are identified, the entire event is kinematically fitted. We 
reconstruct the rj candidates from their decays to 77, 7r + 7r~7r°, and 3tt°. We did not employ 
the decay mode r] — > 7r + 7r~7 because of its small branching fraction (£> = [4.69 ±0.10]% [5j) 
and large backgrounds, primarily from T(nS) — > 7r + ir~T(mS). 

Photon candidates must be detected in the central region of the calorimeter (| cos#| < 
0.81), must not be aligned with the initial momentum of a track, and should have a lateral 
shower profile consistent with that of a photon. Neutral pion candidates (except in the decay 
f] — > 37T°, where we only look for six photon candidates) are reconstructed from a pair of 7 
candidates required to have 77 invariant mass between 120 and 150 MeV. 

Monte Carlo (MC) samples were generated for generic Y(2S, 3S) decays using the routine 
QQ and for T(nS) — > (r]/n )T(mS) and dipion transitions between T states using the 
package EvtGen |l2j. The final T(mS) state was taken to decay to e + e~ or A 
GEANT-based 13( detector simulation was used. These samples, as well as off-resonance 
T(2S) data, are useful both for validating background suppression methods and as possible 
background sources. In calculating branching fractions from data, we take B[T(1S) — > 
e+e-) = B[?(1S) -> n+yT] = 0.0248 ± 0.0005 [5] and B[T{2S) -> e+e"] = B[T(2S) -> 
= 0.0203 ± 0.0009 [14| based on the more accurately measured yU + /i~ branching 
fractions and assuming lepton universality. 

The T(nS) — ► (i]/ir )T(mS) MC samples were generated with rj and 7r° decaying through 
all known decay modes. These decays proceed via a P-wave, and hence are described 
by a matrix element (e^x ej) ■ P( v /n°) in the nonrelativistic limit (here * denotes complex 
conjugation), with e/j the polarization vectors of the final and initial T. The 9 distribution 
for the final-state leptons in T(mS) — > £ + £~ then is 1 — (1/3) cos 2 9, and was used in all signal 
MC samples for T(nS) — > (r]/ii )T(mS). For T(nS) — > wT (mS) it was assumed that the 
T(mS) retains the polarization of the initial T(nS), so the lepton angular distribution for 
T(mS) -^£ + r is l + cos 2 #. 

As a cross-check, data were analyzed for the known transitions T(nS) —> 7nrT(lS), and 
branching fractions were found in sufficiently good agreement with world averages [5j. We 
looked for systematic differences between detection of T(IS') — > e + e~ and T(IS') — > 
Efficiencies for the two modes can differ as a result of the requirement on cos# e + mentioned 
above. The branching fractions calculated from T(1S) — > e + e~ and T(1S) —* [i + fi~ were 
found to be equal within statistical uncertainty, and consistent with those obtained from 
recoil mass spectra without requiring final leptons. 

Kinematic fitting was used to study the decays T(nS) — > (rj / tt°)T (mS) . The two tracks 
selected as leptons, including photon bremsstrahlung candidates within 100 mrad of the 
initial lepton direction, were constrained to have the known masses of T(mS) with a resultant 
reduced \ 2 {x 2 /d-o.f.), Xr = xj+£- m rec L u rred to be less than 10. (For off-resonance data 
the dilepton masses were reduced by an amount equal to the initial M[T(nS)] minus the off- 
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resonance center-of-mass energy.) The sum of the four-momenta of these two fitted tracks, 
including photon bremsstrahlung candidates as well as the decay products of the 77/71-°, 
were further constrained to the initial T(nS) four-momentum, with a reduced x\ = xlvTm 
required to be less than 10, or 3 for (77/71-°) — > 77 to help suppress doubly radiative Bhabha 
events. Some of these Bhabha events can give small fitted Xevt m i but have photon momenta 
shifted by relatively large amounts compared to signal events. To further suppress such 
events, two-photon "pull" masses, defined as (fitted-measured)/<7, where a is the two-photon 
mass resolution, were chosen on the basis of signal MC and off-resonance data (containing 
the doubly radiative Bhabha contribution) to lie between -2 and 3. Over 99% of the signal 
MC events for all transitions satisfy this criterion. All particles were also required to have 
common vertices in the above two constrained fits, with reduced X1+1- v < 30 required for 
the dilepton vertex, and reduced x!vt v < 30 required for the full event vertex. 

For T(2S) — ► T]T(1S), the photons from rj — > 77 have energies E 1 = (281 — 64 cos 9*) 
MeV, where 9* is the angle between the photon in the 77 center-of-mass and the 77 boost, so 
217 < Ej < 345 MeV. Choosing 200 < E 1 < 360 MeV then eliminates background from 
T(2S) — > JXbJ —> 77^(15') with little effect on the 77 — > 77 signal. Using the T(2S) — > 
tjT(1S) MC sample, the 77 candidate mass distribution was fitted to the sum of a double 
Gaussian and a linear background. Constant background gave a worse fit because of the 
kinematic limit at M[T(2S)} — M[T(1S')] = 563 MeV. The fitting range was chosen to be 
533 to 563 MeV: roughly symmetric about the 77 peak {M{rj) = 547.51 ±0.18 MeV with 
upper boundary at M[T(2S)] — M[T(LS)] above which few events are expected or observed. 
The difference between fits with linear and flat backgrounds was found to be insignificant 
compared with the systematic error associated with fitting range. The double Gaussian 
parameters included a narrow width o\ = 0.9 MeV, a wide width 02 = 2.1 MeV, area of 
second peak 20% of total, and mean of second peak 0.14 MeV below the first. 

The mass distribution for the sum of the 77 modes 77, 7r + 7r~7r°, and 37r° in data (up- 
per plot, Fig. [TJ) shows a clear peak near M(rj). We fit data points to the sum of the 
double Gaussian with floating area but fixed shape obtained from signal MC and a lin- 
ear background. The 7r + 7r~7r° and 3tt° decay modes each contribute two events near the 
peak and none elsewhere. The combined fitted peak corresponds to a branching fraction 
B[T(2S) -»■ rjT{LS)] = (2.1±g;2) x 10~ 4 . Defining the significance N a as V-2Alog£, where 
£ is the likelihood, the difference between fits with and without signal yields a statistical 
significance of 5.3 standard deviations. 

In searching for Y(3>S) — » 77(— > 77)T(15') transitions, we suppress backgrounds from 
cascades involving intermediate Xb(l-P>2P) states by requiring one photon to have 500 < 
E 1 < 725 MeV and the other to have 140 < E 2 < 380 MeV. Signal photons satisfy £ 7 = 
(435 — 350 cos#*) MeV, so about 2/3 of them are retained by these choices. Small differences 
with respect to T(2S) — > r]T(lS) include (a) an 77 fit range 523-573 MeV and (b) a flat 
background, found here to be sufficient to describe MC and data. The best fit to signal MC 
shape and the 90% confidence level (CL) upper limit are shown in the lower plot of Fig. [TJ 
(No events were observed in the regions included in the fit but not shown in Fig. [TJ) 

For T(2S) —>■ 7r°T(15), the photons from 7r° — > 77 have energies _E 7 = (274 — 266 cos#*) 
MeV, so 8 < £ 7 < 540 MeV. The choice 200 < £ 7 < 360 MeV for both photons, made 
to eliminate background from T(25) — > •"fXbJ 77Y(1S), then retains about 30% of the 
7T° — > 77 signal. A fit of the M(77) distribution in the data (using the signal MC double- 
Gaussian shape and uniform background) is shown in the top plot of Fig. [2J Details of 
this and other limits, as well as of the T(25) — > r]T(lS) signal, are shown in Table [TJ 



5 



1850508-001 



II 

J — h- 


-I 1 1 1 

Y(2S)^r,Y( 
- , rrt 


,S '/1 


1 — 1 — 1 

\ 


ii on 


- 


1) 1 


-H 1 h 

(1 

( 


T 

» 


I 1 "" 

(3S)- 

< 


■ 


j — i — 
1S) 


H H— H 

» 


» ■ 



533 543 553 563 

Invariant mass of r\ candidate (MeV) 



FIG. 1: Events per MeV vs. invariant mass of candidates for T(25) — ► nT(lS) (top) and T(35) — > 
i]T(lS) (bottom). The sum of the modes n — > 77, n — > tt + tt~tt°, and n — > 37r° is shown. In the 
top figure the solid curve corresponds to the total fit, involving a signal of I3.9I3J even ts above 
background (dashed line). In the bottom figure the solid curve corresponds to a best fit with signal 
MC shape, while the dotted curve corresponds to a 90% confidence level (CL) upper limit. 



For all 7T° transitions, MC simulations indicate a constant function is adequate to describe 
the background. Efficiency differences between decay modes are typically due to details of 
photon acceptance. 

For T(35') — > tt°T(1S), where signal photons from tt° — > 77 satisfy E 1 = (429 — 
385 cos#*) MeV, the same ranges of (Ei,E 2 ) are chosen as for T(35*) — > 7]T(1S). For 
T(3S') — > ir T(2S), we suppress backgrounds from cascades involving intermediate x?>(2-P) 
states by excluding photons with 60 < E 2 < 130 MeV and 190 < E 1 < 260 MeV. Here, the 
signal photons satisfy E 1 = (164— 149 cos 9*) MeV, so about 40% are retained. No signal is 
seen in any of these 7r° transitions (Fig. [2]). 

Systematic errors are shown in Table [III Other contributions investigated and found to 
be negligible were (i) cross feeds among rj modes, (ii) signal shape, (iii) background shape, 
(iv) triggering details, and (v) differences in e//x reconstruction. The dominant sources of 
systematic uncertainties are described below. 

1. Bhabha event suppression: Uncertainties for all processes will arise from our Bhabha 
event suppression requirement. Although it is applied only to 77 modes, it will affect 
not only tt° transitions but also those with rj, whose 77 decays dominate our analyses 
statistically. To probe this uncertainty, we consider the separate sample of those events with 
cos# e + > 0.5 which were removed by the Bhabha suppression requirement. The resultant 
£>[Y(2S') — > r)(— > 77)T(1S')] is consistent with our nominal result. Averaging the two gives a 
deviation of 9% which we take as a possible systematic uncertainty due to this requirement. 
We then propagate this estimated uncertainty to the rest of the decay modes with suitable 
weight for the fraction of the decay due to 77. 

2. Kinematic fitting: To probe any systematic bias introduced by our kinematic fitting 
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FIG. 2: Best fits to two-photon invariant mass distributions with signal MC shapes (solid curves 
are the results of total fits) and 90% CL upper limits (dotted curves) for T(2S) — > n°T(lS) (top), 
T(35) -» tt°T(15) (middle), and T(35) -» vr°T(25) (bottom). 



TABLE I: Efficiencies, events in data, and product branching fractions B x Be, where B = 
B[T(nS) -> (r//7r°)T(m5)], and = B[T(15 r ) -► £ + £~] = 4.96% or B[T(25) -> = 
4.06% (^ + £ _ = e + e~ + /j + /j~). Efficiencies are based on MC samples generated with standard 
77 and 7T° branching fractions and with B[T(mS) — > e + e~] = B[T(mS) — > = 50%. Decays 

involving n are based on combined 77, 7r + 7r~7r°, and 3tt° modes. 
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procedure, we look at events with very similar topology to our signals: T(2S) — > 7% 6 (lPj), 
Xft(l-Pj) — > 7T(1S'), T(1S) — * where J = 1 or 2. We use the same analysis requirements 
as for i] — > 77 but relax the requirements on i? 7 in order to accept two-photon cascades 
through states. Varying the requirement on xlvTm from none to xlvTm < 3, we observe 
a maximum deviation of 7% in this product of branching fractions which we assign as a 
possible source of systematic uncertainty. 

3. rj/n reconstruction: We assign 5% per tt° or 77 decaying into two photons based on 
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4. Fit ranges: Uncertainties due to fit ranges differ for different final states. To estimate 
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TABLE II: Systematic errors, in percent, on branching fractions for T(nS) — ► (a) rjT(lS); (b) 
7r°T(l<S); (c) tt°T(2S). All errors are assigned symmetrically. Decays involving r\ are based on 
combined 77, 7r + 7r~7r , and 3ir° modes. The last line (d) includes systematic errors. 
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them, we prepare many MC samples in which points are randomly scattered around best-fit 
values from data (signal plus background), bin- by-bin according to a Poisson distribution. 
We then fit them with the fit range boundaries symmetrically changed by ±5 MeV for 
T(2S) — > T]T(1S). In T(3S') — > r]T(lS) as well as in T(nS) — * 7r°T(m5'), where wider 
kinematic ranges are available, the fit range boundaries are symmetrically changed by ±10 
MeV. We assign variations of averages of these fitted yields as possible systematic shifts. 
Combining the effects from the systematic errors linearly with the results already listed, we 
find the results shown in the last line of Table [Til 

To summarize, we have observed for the first time a process involving 6-quark spin-flip, 
with B[T(2S) -> r/r(l£)] = (2.1^ ± 0.3) x 10~ 4 . The statistical significance of the signal 
is 5.3 o. The result is about 1/4 of the value one would predict on the basis of Eq. ([I]), 
indicating either a shortcoming in the description of two-gluon hadronization into an 77 or 
a fundamental suppression of the chromomagnetic moment of the b quark. In addition, 
we have set 90% CL upper limits on other pseudoscalar transitions summarized on the 
bottom line of Table [TTJ The limit on B[T(3S) — ► rjT(lS)] is about a factor of two below 
that predicted from Eq. ([I]), while the limits on the transitions T(2S,3S) — > n°T(lS) are 
consistent with the estimates of Eq. (j2J) . 
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